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ABSTRACT 

Comparison of quasar absorption spectra with laboratory spectra allow us to probe 
possible variations in the fundamental constants over cosmological time-scales. In a 
companion paper we present an analysis of Keck/HIRES spectra and report possible 
evidence suggesting that the fine structure constant, a, may have been smaller in the 
past: Aa/a = (-0.72 ± 0.18) x 10~ 5 over the redshift range 0.5 < z < 3.5. In this 
paper we describe a comprehensive investigation into possible systematic effects. Most 
of these do not significantly influence our results. When we correct for those which do 
produce a significant systematic effect in the data, the deviation of Aa/a from zero 
becomes more significant. We are lead increasingly to the interpretation that a was 
slightly smaller in the past. 

Key words: atomic data - line: profiles - instrumentation: spectrographs - methods: 
data analysis - techniques: spectroscopic - quasars: absorption lines 



1 INTRODUCTION 

The idea of varying fundamental constants is not a new 
one: Milne (1935, 1937) and Dirac (1937) independently sug- 
gested a varying Newton gravitational constant. More mod- 
ern theories, such as superstring theory, naturally predict 
variations in fundamental coupling constants (e.g. Forgacs 
& Horvath 1979; Marciano 1984; Barrow 1987; Li & Gott 
1998), strongly motivating an experimental search. Sev- 
eral authors have demonstrated the idea (Varshalovich & 
Potekhin 1995 and references therein), first used by Bahcall, 
Sargent & Schmidt (1967), of using spectra of high redshift 
gas clouds seen in absorption against background quasars 
(QSOs) to constrain possible variation of the fine struc- 

2 

ture constant, a = j—, over cosmological time-scales. This 
method is based on a comparison of the observed transition 
wavelengths of alkali doublets and recent applications of it 
have yielded constraints on Aa/a = (q 2 -qo)/qo ~ 5x 10~ J 
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nia and the National Aeronautics and Space Administration. The 
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where a z and ao are the values of a at the absorption cloud 
redshift z and in the laboratory respectively (e.g. Cowie & 
Songaila 1995; Ivanchik, Potekhin & Varshalovich 1999; Var- 
shalovich, Potekhin & Ivanchik 2000). Recently, we have in- 
creased the level of precision obtained with the alkali doublet 
method to Aa/a — (— 0.5± 1.3) x 10 -5 using higher quality 
data (Murphy et al. 2000c). 

However, a new method offering an order of magnitude 
improvement in precision was suggested by Dzuba, Flam- 
baum & Webb (1999a,b) and was demonstrated by Webb et 
al. (1999, hereafter W99). The increase in precision derives 
from the fact that the relativistic corrections to the energy 
levels of different ions vary from species to species by up 
to an order of magnitude. A comparison of the observed 
Mgl, Mgn and Fell transition wavelengths in 30 absorption 
systems allowed W99 to tentatively suggest that the mean 
value of Aa/a = (-1.09 ± 0.36) x 10" 5 over the redshift 
range 0.5 < z < 1.6: the fine structure constant may have 
been smaller in the past. 

The present work is the companion paper of Murphy 
et al. (2001a, hereafter MOla) in which we discuss our tech- 
niques in detail and present the results of recent work. We 
summarise the main points here. We obtained Keck/HIRES 
spectra of 28 QSOs and analysed 49 absorption systems ly- 
ing along their lines of sight. These systems included those 
Mg n/Fe II systems considered in W99 and also a higher red- 
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shift sample of damped Lyman-a systems (DLAs) contain- 
ing absorption lines of ions such as Mgl, Mgn, Alii, Aim, 
Sin, Cm, Fen, Nin and Znn, so that our redshift range 
now covers 0.5 < z < 3.5, corresponding to look-back times 
from 5 to 12 billion years (Ho = 68 kms - Mpc - , Qm = 0.3, 
SI a = 0.7). Our absorption systems fell conveniently into 
two sub-samples: the Mgll/Fell systems comprised the low 
redshift sample (z = 1.0) and the DLAs comprised the high 
redshift sample (z = 2.1). There were, however, 2 Mgll/Fen 
absorbers within the high-z sample that lay at low redshift 
but we have included them in the high-z sample since the ob- 
servations and data reduction for the two sub-samples were 
carried out independently by two different groups. 

The analysis is based on the following equation for the 
wavenumber, u z , of a given transition at a redshift z; 

uj z = ljo + qix + q2y , (1) 

where uio is the laboratory rest wavenumber of the transi- 
tion, qi and Q2 are relativistic coefficients, x = (a z /ao) 2 — 1 
and y = (a z /ao) 4 — 1. Note that the second and third terms 
on the right hand side only contribute if Aa/a 7^ 0. The 
values of q\ and §2 have been calculated using many-body 
techniques in Dzuba et al. (1999a, b, 2001) and we tabulate 
those values relevant to our analysis in table 1 of MOla. 

The values of u>o must be known to a high precision 
in order to place stringent constraints on Aa/a. A preci- 
sion of Aa/a ~ 10~ 5 corresponds to an uncertainty on wo 
~ 0.02 cm -1 (~ 0.3 kms _1 )|^]. This corresponds fairly well 
to the data quality since we can centroid an absorption fea- 
ture to approximately one tenth of the HIRES resolution, 
~ 7 kms - . Such precision has not been previously achieved 
for most resonance transitions of ionized species and so mea- 
surements of u>o for all relevant species were carried out by 
several groups (Nave et al. 1991; Pickering, Thorne & Webb 
1998; Pickering et al. 2000; Griesmann & Kling 2000) using 
laboratory Fourier transform spectrometers. We have sum- 
marised and tabulated these measurements in MOla. 

Knowing all values of loq, q\ and 52, we then decon- 
volve each absorption system into individual Voigt compo- 
nents and perform a non-linear least squares \ 2 minimiza- 
tion to find the best-fitting value of Aa/a at each absorp- 
tion redshift. For the low-z sample we confirm the W99 re- 
sult, our new results giving a weighted mean of Aa/a = 
(-0.70 ±0.23) x 10~ 5 . For the new, high-z sample, we find a 
similar result: Aa/a = (-0.76±0.28) x 10~ 5 . Both results in- 
dependently suggest a smaller value of a in the past. Taking 
the sample as whole we find Aa/a = (-0.72 ±0.18) x 10" 5 , 
a 4cr result. These new results are also summarised in Webb 
et al. (2001, hereafter W01). 

Given the potential importance of such a result for fun- 
damental physics, a thorough investigation of possible sys- 
tematic problems in the data or in the analysis is essential, 



t We can see this as follows. We can ignore the term in 92 since 
these coefficients are typically an order of magnitude smaller than 
the q\ coefficients. Consider two transitions, one with a low atomic 
mass, the other a high atomic mass. Typically, the difference be- 
tween the <ji coefficients for a light species (e.g. Mg 11) and a heavy 
species (e.g. Fell) is ~ 1000 cm± Equation [l| applied to these 
two transitions, then shows that a precision of ojo of ~ 0.02 cm -1 
is required in order to detect a variation of Aa/a ~ 10 . 



and this is the purpose of the present paper. In Section 2 
we list all effects that one might consider as being able to 
mimic a systematically non-zero Aa/a. We are able to ex- 
clude many of them with general arguments. Sections 3-7 
are then given over to a more detailed analysis of five specific 
effects: wavelength mis-calibration, ionic line blending, at- 
mospheric dispersion effects, differential isotopic saturation 
and possible isotopic variation. Section 8 describes another 
test which is sensitive to a simple monotonic wavelength 
distortion, specifically for the high-z sample. We make our 
conclusions in Section 9. 



2 SUMMARY OF POTENTIAL SYSTEMATIC 
ERRORS 

In MOla and W01 we found that values of Aa/a were 
systematically negative and that this trend was evident in 
two different samples of data occupying different redshift 
regimes. We note that the low-z sample is more suscepti- 
ble to systematic errors than the high-z sample. The Mgn 
lines act as anchors against which the larger Fe n shifts can 
be measured. All of the Fe 11 lines have similar magnitude 
shifts in the same direction (i.e. qi is large and positive in all 
cases). Furthermore, all the Fen lines lie to the blue of the 
Mgn doublet. If some systematic error in the wavelength 
scale is present, which effectively compresses the spectrum, 
this would mimic a negative shift in the mean value of Aa/a 
(if it were not degenerate with redshift). The high-z sample 
is far more complex in this regard since it contains ions not 
only with positive coefficients, but also with negative qi co- 
efficients (i.e. all Crn lines and Nin A1741 and A1751). Also, 
the shifts for the various ions are of various magnitudes mak- 
ing it more complicated to predict the effect of a systematic 
error on our measured values of Aa/a. 

2.1 Laboratory wavelength errors 

Errors in the values of u>o will lead directly to errors in any 
single determination of Aa/a. If, for example, the labora- 
tory wavelengths of the Mgn lines were all shifted to the 
blue or red then this would systematically bias the values 
of Aa/a for the low-z sample. In MOla we summarise the 
measurements of all values of ujo- Most measurements were 
repeated independently and good agreement was obtained 
in those cases. A typical measurement accuracy from these 
laboratory measurements ujo of is ~ 2 x 10~ 3 cm -1 , which 
could only introduce a maximum error in Aa/a around an 
order of magnitude below the observed deviation from zero. 
In practice, since we use several transitions (particularly for 
the high-z sample), any associated errors are likely to be 
smaller. We therefore consider it unlikely that the system- 
atic shift in Aa/a is due to errors in the values of u>o. 

As also noted in MOla, some laboratory spectra were 
calibrated with the Aril lines of Norlen (1973) while others 
were calibrated with those of Whaling et al. (1995). The 
Whaling wavenumbers are systematically larger than those 
of Norlen such that Auj = 7 x 10~ 8 tj. To deal with this 
discrepancy we have re-scaled those wavenumbers which had 
been calibrated using the Whaling scale (Aln, Aim, Sin) 
to the Norlen scale. We stress that the difference between 
the Norlen and Whaling wavenumbers is a scaling only, and 
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not a distortion. Therefore, is it not possible to introduce 
an apparent deviation in Aa/a from zero, since small zero- 
point offsets in the wavelength scale are degenerate with 
redshift. 



expect to find Aa/a ~ 0.6 x 10 -5 in the absence of any real 
variation in a. 

We consider this possibility in detail in Section M 



2.2 Heliocentric velocity variation 

During a one hour exposure, representative of the quasar 
integrations for the sample we have used, the heliocentric 
velocity may change by as much as ~ 0.1 kms -1 . This will 
act to smear any spectral feature: the instantaneous spec- 
trum is convolved with a top hat function in velocity space. 
However, since we fit a redshift parameter to the spectrum 
when determining Aa/a then a velocity space smearing is 
completely absorbed into the redshift parameter. Heliocen- 
tric smearing will have no effect on measurements of Aa/a. 



2.3 Differential isotopic saturation 

The ions used in our analysis, with the exception of those 
of Al, have naturally occurring isotopes. Thus, each absorp- 
tion line will be a composite of absorption lines from all the 
isotopes of a given species. In all cases other than Mg and 
Si, the laboratory wavenumbers we have used correspond 
to composite values. The composite wavenumbers will only 
strictly be applicable in the optically thin regime (linear part 
of the curve of growth). As the column density increases, the 
strongest isotopic component begins to saturate and the line 
centroi d wi ll shift according to the natural abundances (cf. 
Section 2.4) of the other isotopes. 

The mass isotopic shift, being the dominant splitting 
mechanism for light atoms, is proportional to u>o/m 2 for 
m the atomic mass. The wavelength separation between 



the isotopes is therefore largest for the lighter (low qi) 
ions. Isotopic separations have been measured for the Mg I 
A2853 transition by Hallstadius (1979), for Mgn A2796 by 
Drullinger et al. (1980) and the isotopic structure of all Mg 
lines is presented in Pickering et al. (1998). However, similar 
measurements do not exist for the other transitions of inter- 
est. For Si we used estimates of the isotopic wavenumbers 
based on a scaling from the basic Mg structure by the mass 
shift. For all other species, measured (composite) wavenum- 
bers were used. See MOla for full details. 

We discuss and quantify this effect in Section 



2.4 Isotopic abundance variation 

In the above we have assumed that the isotopic abundances 
in the absorption systems are equal to terrestrial values (Ta- 
ble ^|). However, if the isotopic abundances at high redshift 
are significantly different then this may lead to apparent 
shifts in the absorption line centroids and this would cause 
a systematic error in Aa/a. Again, the most susceptible 
transitions are those of Mg and Si since (i) these ions will 
have the largest isotopic separations (i.e. m is low) and (ii) 
transitions in these atoms are prominent in our data. If, 
for the low-z sample, the 25 Mg and 26 Mg abundances were 
zero, we would incorrectly infer a more positive Aa/a by 
assuming terrestrial isotopic abundances. In the simple case 
of finding Aa/a from Mgn A2796 and FenA2344, we would 



2.5 Hyperfine structure effects 

Hyperfine splitting of energy levels occurs in species with 
odd proton or neutron numbers. Different hyperfine com- 
ponents will have different transition probabilities but the 
composite wavelength of a line will be unchanged by the 
splitting (i.e. the centre of gravity of the hyperfine compo- 
nents is constant with increased splitting). However, a sim- 
ilar differential saturation effect wil l oc cur for the hyperfine 
components as discussed in Section 2.3 for the isotopic com- 
ponents. 

The worst case here is that of 27 A1 - the only stable 
isotope of Al - since it has an odd proton number and a 
very large magnetic moment. In the laboratory experiments 
of Griesmann & Kling (2000), the hyperfine structure of 
the Aim doublet was clearly resolved and can be seen to 
be quite prominent. However, as noted in MOla, we only 
include Aim lines in the analysis of 6 absorption systems 
of our high-z sample and in these cases we see very little or 
no saturation. The effect of removing the individual Aim 
transitions from those systems on Aa/a can be seen in Fig. 
4. In the case of Alii A1670, both the ground and excited 
state have zero spin and the p-wave hyperfine splitting in 
the excited state is very small (cf. 25 Mg in Pickering et al. 
1998). Thus, we do not expect any Al hyperfine saturation 
effects on Aa/a. 

All other dominant isotopes in our analysis have even 
proton numbers and where odd isotopes occur these have 
very low abundances (see Table ^). The 25 Mg hyperfine 
structure has already been taken into account in our analy- 
sis but the hyperfine separations are so small that neglecting 
it would make no difference to Aa/a even upon saturation. 
The odd isotopes of Si, Cr, Fe, Zn and Ni have very low 
abundances. Also, Si, Cr, Fe and Ni will have very small hy- 
perfine splitting due either to low magnetic moments or to 
the nature of the ground and excited state wavefunctions. 

We therefore expect that differential hyperfine struc- 
ture saturation has not effected Aa/a in our analysis. An 
unlikely interpretation, given the above, is cosmological evo- 
lution of the hyperfine structure constants but very large 
relative variations would be required. 

Another possibility is that the populations of the two 
hyperfine levels of Aim are not equal in the absorption 
clouds. Since the gas clouds have very low density, the equi- 
librium between the two hyperfine levels will be maintained 
predominantly by interaction with cosmic microwave back- 
ground (CMB) photons. For the AlmA1862 transition, this 
implies a lower limit on the relative populations of the two 
levels of exp( — A^/A;bTcmb) ~ 0.9, leading to a shift in the 
line centroid of - 0.01cm" 1 (~5x 10" 4 A). This is a po- 
tential problem for the 6 systems in which we fit Aim lines. 
We do not consider this effect in detail but we note that 



Fig. 4 in Section 4.2 shows that potential lack of thermal 



equilibrium has no significant effect on Aa/a. 
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2.6 Magnetic fields 

Shifting of ionic energy levels due to large magnetic fields 
may also be a possibility. If very large scale magnetic fields 
exist in the intergalactic medium and quasar light is polar- 
ized, then this could result in correlated values of Aa/a in 
neighbouring regions of the universe. However, the QSOs 
studied in MOla and W01 cover a large region of the sky. 
Furthermore, even in Abell clusters, intra-cluster magnetic 
field strengths are typically ~ /j,G (e.g. Feretti et al. 1999), 
roughly 9 orders of magnitude below the strength required 
to cause substantial effects. We consider this possibility to 
be very unlikely. 



2.7 Kinematic effects 

An assumption in the analysis of MOla and W01 is that the 
absorption velocity structure is the same in different species. 
Of course, if there were some non-zero velocity between, say, 
the Mgn and Fell components of a given cloud then we 
would incorrectly determine Aa/a. However, it is difficult 
to imagine a mechanism by which the Mg II is systematically 
blue or redshifted with respect to the Fen when averaged 
over a large sample of absorption systems. Suppose for ex- 
ample we have an asymmetric, expanding (or contracting), 
rotating cloud in which all the Fe II and Mg II lie at different 
'radii'. Any large sample of sight lines through such a cloud 
would yield an average value of Aa/a — 0. Furthermore, 
if this was a major effect then the observed scatter in the 
Aa/a points would be greater than expected on the basis 
of our Aa/a error estimates, and in MOla we showed this 
not to be the case. No matter how contrived the velocity 
structure and ionic segregation, when averaging over a large 
number of values of Aa/a, it is not feasible to generate a 
systematic offset from zero. 



2.8 Wavelength mis-calibration 

There may be inaccuracies in the wavelength calibration 
software used to reduce the spectra. Also, human error can- 
not be ruled out. W99 stated that a drift in the wavelength 
calibration drift by an amount corresponding to roughly 
~ 2.5 times the mean wavelength-pixel residuals over the 
range of the CCD could result in a significant error in the 
law-z results. A considerably more complicated form of the 
mis-calibration would be required to produce a systematic 
error in the high-z results. 

The wavelength calibration of the QSO CCD images 
is done by comparison with Thorium-Argon (ThAr) lamp 
exposures taken before and after the QSO frames. We con- 
sider errors in the laboratory wavelengths of the selected 
ThAr lines to be negligible as their relative values are known 
to much greater precision than our values of ljq (Palmer & 
Engleman Jr. 1983). However, ThAr line mis-identifications 
may lead to serious mis-calibration of the wavelength scale 
over large wavelength regions. Such mis-identifications can 
also be applied to the rest of the spectra taken over the 
same series of observations (i.e. applied to the spectra of 
many QSOs). Thus, a priori, we cannot rule out the possi- 
bility that the wavelength scale has been set improperly in 
this process and that a systematic shift in a has not been 
mimicked. Such a potential effect needs careful investigation. 



In Section |§] we discuss and quantify this effect for our 
data set. 



2.9 Air— vacuum wavelength conversion 

Most ThAr line lists are presented as air wavelengths. The 
usual data reduction procedure is to carry out the wave- 
length calibration fits to the data using air wavelengths 
and then convert the calibrated spectrum to vacuum wave- 
lengths. For example, IRAF^] uses air wavelengths for a se- 
lection of lines from the Thorium spectral atlas of Palmer 
& Engleman Jr. (1983) and from the Argon lines of Norlen 
(1973). The same list is used in makee - the HIRES data 
reduction package written by T. Barlow, makee was used 
to calibrate 12 of our high-z sample and the remainder were 
calibrated within iraf. The latter were actually calibrated 
with a set of vacuum wavelengths, A vac , that iraf states is 
derived from the air wavelengths, A a i r , by using the Edlen 
(1966) formula for the refractive index, n, 



10 s [n(A vac ) - 1] = 8342.13 + 



2406030 



+ 



15997 



(130 -a 2 ) (38.9 -a 2 ) 



(2) 



at 15°C and atmospheric pressure. Instead of using a = 
10 4 /A V ac, IRAF makes the approximation, a « 10 4 /A a ir- ma- 
kee converts the final air wavelength calibrated QSO frames 
to vacuum using the Cauchy dispersion formula (Weast 
1979), 

0.3555 



1 O OQQ 

10 7 [n(A air ) - 1] = 2726.43 + 1Q + 



10- 16 A 4 . 



(3) 



at 15° C and atmospheric pressure. A vac and A a i r are both 
measured in Angstroms in the above equations. This differ- 
ence begs the question, what is the absolute and relative 
accuracy of these conversion formula? 

First let us refine this question. The measurements re- 
ported by Palmer & Engleman Jr. and Norlen were all car- 
ried out in vacuum. The air wavelengths they quote are cal- 
culated using the Edlen formula. Thus, iraf and makee will 
only produce 'correct' vacuum wavelength scales if they in- 
vert the Edlen formula to calculate the refractive index. As 
stated above, this is not the case. Therefore, to re-state the 
above question, what error does the iraf and makee air- 
vacuum wavelength conversion introduce into the resulting 
QSO wavelength scale? To illustrate the answer to this we 
plot the various conversion formulae commonly used in the 
literature, including those discussed above, in Fig. 1. We 
have also included the orig inal Edlen (1953) formula^] for 
comparative purposes. Peck & Reeder (1972) published a 
fit to revised refractive index data, proposing a 4-parameter 
formula: 

5791817 167909 



10 s [n(A vac ) - 1] = 



+ ■ 



(238.0185 -a 2 ) (57.362 - a 2 ) 



(4) 



S iraf is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ The Edlen (1953) formula is of the same form as the 1966 
formula (equation Eh but has slightly different parameters: 6432.8, 
2949810, 146, 25540 and 41 in order from left to right, numerator 
to denominator. 
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Figure 1. Comparison of the various air dispersion formulas used 
in the literature (all for 15° C and atmospheric pressure). We use 
the Edlen (1966) formula as our reference since the atlases of 
Palmer & Engleman Jr. (1983) and Norlen (1973) use this to 
convert their vacuum wavelengths to air wavelengths. The top 
panel shows the dispersion for the Edlen (1966) and Edlen (1953) 
formulas together with that of Peck & Rccdcr (1972) and the 
Cauchy formula (Weast 1979). We also plot the approximation 
made within I RAF where the Edlen formula has not been inverted 
to convert air wavelengths to vacuum (denoted "Edlen n.i."). The 
middle and lower panels show the difference between the Edlen 
(1966) refractive index, n^, and that derived from the other for- 
mulae, tlx: An = riE — "x- 

at 15° C and atmospheric pressure with a = 10 4 /A vac (again, 
A vac is in Angstroms). 

From Fig. 1 we see that the Cauchy formula seriously 
deviates from the other formulas. Furthermore, the differ- 
ence between the Edlen and Cauchy dispersions is strongly 
wavelength dependent. Thus, spectra calibrated with MA- 
kee may have systematically non-zero values of Aa/a due 
to this distortion. As we stated in MOla, we have corrected 
the 12 effected spectra in our (high-z) sample by converting 
the wavelength scale back to air wavelengths using the (in- 
verted) Cauchy formula and then back to vacuum using the 
(inverted) Edlen formula. From Fig. 1 we can also see that 
there is only a small distortion to the wavelength scale in- 
troduced due to the approximation made within iraf. This 
difference is below our level of precision. 



2.10 Temperature changes during observations 

The refractive index of air within the spectrograph depends 
on temperature (and also on pressure but this is a smaller 
effect). If a QSO spectrum is calibrated with only a single 
ThAr exposure taken at a different spectrograph temper- 
ature, we expect a systematic mis-calibration of the QSO 
frame's wavelength scale. Consider a QSO frame taken at 
15° C which is calibrated with a ThAr exposure taken at 0°C 



and similar pressure. Any wavelength separation in the QSO 
exposure will be overestimated. For example, the difference 
between the wavelength separations of vacuum A4000 and 
A7000 A at 15°C and 0°C is 43 x 1CT 3 A (using the tables 
in Weast 1979). A comparison with values of A A in table 1 
of MOla shows that this would result in a significantly pos- 
itive Aa/a for the low-z sample, the situation not being so 
clear for the high-z sample. 

Note that the above effect can only mimic a system- 
atically non-zero Aa/a if the spectrograph temperature is 
systematically higher or lower for the ThAr frames com- 
pared with the QSO frames. This is possible if the ThAr 
exposures were always taken before or after the QSO frames 
and the temperature evolves monotonically throughout the 
night. However, the effect is greatly reduced if ThAr expo- 
sures are taken near the time of the QSO observations. This 
was the case in our observations. We have used image header 
information to calculate the QSO-ThAr temperature differ- 
ence, AT = (Tqso) - (T Th Ar) , inside HIRES for both the low 
and high-z samples. Here the average of Tqso and TrhAr is 
taken over all exposures for each object. We find mean val- 
ues of AT = 0.04 ±0.02 Kand AT = 0.2±0.1 K for the low 
and high-z samples respectively. Taking into account that, 
on average, rest-frame separations between transitions are 
300 A then we can see that temperature variations could 
not have mimicked any significant shift in Aa/a in either 
sample. 

2.11 Line blending 

The errors in Aa/a presented in MOla and W01 take into 
account errors from signal-to-noise and spectral resolution 
considerations and the velocity structure of the profile fits. 
The errors are also reduced when more lines are fitted simul- 
taneously. However, we have assumed that we have decon- 
volved each absorption system into the correct number of 
velocity components. There may have been weak, interlop- 
ing, unresolved lines which, if the interloping species were 
in the same absorption cloud, could have produced a shift 
in the fitted line wavelengths of all velocity components of 
one or more transitions. 

We distinguish between random blends and systematic 
blends. Random blends may occur if many absorption clouds 
at different redshifts intersect the line of sight to a single 
QSO. A systematic blend will occur when two species are in 
the same cloud and have absorption lines with similar rest- 
wavelengths. Such an effect could mimic a systematic shift 
in a. 

The importance of such an effect diminishes as the num- 
ber of transitions used in each fit is increased. Therefore, 
the effect, if present, would be expected to be smaller in our 
high-z sample. 

In Section |^ we describe the results of a detailed search 
for atomic (and not molecular) interlopers for all transitions 
of interest, and place constraints on their strengths and po- 
sitions. 

2.12 Atmospheric dispersion effects 

As noted in MOla, all of the low-z sample and many of the 
high redshift sample of QSOs were observed before 1996 Au- 
gust at which time an image rotator was fitted to HIRES. 
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Objects observed prior to August 1996 were therefore not 
observed with the slit length (spatial direction) perpendicu- 
lar to the horizon. Atmospheric dispersion leads to a stretch- 
ing of the target spectrum relative to the calibration spec- 
trum. 

An additional consequence is the potential wavelength- 
dependent truncation of the QSO seeing profile on the slit 
jaw edges. This can introduce a wavelength dependent asym- 
metry in the point-spread function (PSF). 

Together, these effects can conspire to mimic a non- 
zero Aa/a. We provide detailed calculations for these effects 
in Section 5 and show that they will tend to produce an 
apparent positive Aa/a, and so cannot mimic our results. 

2.13 Instrumental profile variations 

If the instrumental profile (IP) of HIRES shows significant 
intrinsic asymmetries then we expect absorption line cen- 
troids to be incorrectly estimated. If the asymmetry varies 
with wavelength, this could mimic a non-zero Aa/a. Valenti, 
Butler & Marcy (1995) have determined the HIRES IP for 
several positions along a single order and they find that 
the IP is indeed asymmetric and that the asymmetry varies 
slightly along the echelle order. Valenti et al. (1995), do not 
quantify the asymmetry variation across orders although it 
is likely to be comparable to the variation along the orders. 

As a by-product of the detailed investigation into pos- 
sible wavelength calibration errors, we find that this effect 
is negligible (see Section |§]for details). 

We have therefore eliminated, with some reliability, po- 
tential systematic errors due to: laboratory wavelength er- 
rors, heliocentric velocity variation, hyperfine structure ef- 
fects, magnetic fields, kinematic effects, air-vacuum wave- 
length conversion errors and temperature changes during 
observations. 

In the following sections we investigate in considerable 
detail the remaining effects: wavelength mis-calibration (and 
instrumental profile variations), line blending, atmospheric 
dispersion effects, differential isotopic saturation and iso- 
topic abundance variation. We also consider a simple test for 
any simple but unidentified systematic errors for the high- 
z sample in Section H. 



3 DETAILED ANALYSIS: WAVELENGTH 
CALIBRATION ERRORS 

In this section, we investigate the possibility that wavelength 
calibration errors could have lead to an apparent non-zero 
Aa/a. To quantify this directly, we analyse sets of ThAr 
emission lines in the calibration lamp spectra in the same 
way as each set of QSO absorption lines has been anal- 
ysed. If no calibration error is present then we expect that 
(Aa/a)ThAr = for all clouds (i.e. for all z). 

3.1 Method 

For each QSO spectrum, we selected ~15 A sections of ThAr 
spectra at wavelengths corresponding to the observed ab- 
sorption lines in the QSO spectra. Each ThAr ~15A sec- 
tion contains several (typically up to ~ 10) lines. Error ar- 



rays were generated assuming Poisson counting statistics. 
We fitted continua to these ThAr sections, dividing the raw 
spectra by the continua to obtain normalized spectra. 

The main steps in the procedure are then as follows. We 
select several independent sets [jj] of ThAr lines. Each set of 
ThAr lines therefore corresponds to, and has been selected 
in an analogous way to, one QSO absorption system. 

In order to obtain the best estimate of any wavelength 
calibration errors, the average (Aa/a)xhAr is determined 
for each ThAr set (typically 3-7 sets for each QSO ab- 
sorption system). The process of defining sets in this way 
also provides a direct additional estimate of the errors on 
(Aa/a) T hAr. 

To fit the ThAr spectra, a modified version of vpfit 
was used (to fit Gaussian profiles, see MO la), adopting ThAr 
laboratory wavelengths from Palmer & Engleman Jr. (1983) 
for values of wo in equation [j]. The qi and q2 coefficients for 
the QSO absorption lines are applied to the corresponding 
ThAr lines. That is, we treat the set of ThAr lines as if they 
were the QSO lines. The free parameters involved in each 
fit were line width (assumed the same for all ThAr lines 
in a given set, although the results were insensitive to this 
assumption), 'redshift' (numerically close to zero) and peak 
intensity. 

3.2 Results and Discussion 

We applied the above method to each available ThAr spec- 
trum corresponding to each QSO absorption cloud consid- 
ered and the results are illustrated in Fig. 2. Our first ob- 
servation is that the scatter in these results is an order 
of magnitude less than that of the QSO absorption results 
themselves. This clearly demonstrates that wavelength mis- 
calibration of the CCD has not mimicked a shift in a to any 
significant degree. The results above also show that the data 
quality permit a precision of Aa/a ~ 10 -7 . 

Secondly, many points in Fig. 2 deviate significantly 
from zero. There are three possible reasons for this: weak 
blended emission lines, ThAr line mis-identifications, and 
errors in the ThAr laboratory wavelengths. Inspecting the 
ThAr spectra does indeed reveal typically ~20 easily iden- 
tifiable, but very weak, lines per 15 A and we would expect 
this to add to the scatter in the final results for (Aa/a)ThAr- 
Furthermore, this effect would not produce correlated de- 
viations in (A«/a)xhAr, as observed. Consistent with this 
interpretation, we found that the value of x 2 per degree of 
freedom for a fit to a given set of ThAr lines was 2> 1. 
Note that if ThAr line mis-identifications had occurred, one 
may expect systematically correlated deviations. Finally, we 
do not expect any significant errors i n th e ThAr laboratory 
wavelengths (as discussed in Section < 



On inspection of some of the ThAr spectra, we noted 
that some ThAr lines towards the blue end of the spectrum 
were slightly asymmetric. This may indicate a degradation 
in the polynomial fits near the edges of the fitting regions. 
Ther e may also be intrinsic asymmetries in the IP (Section 
|2.13] ). If such asymmetry variations exist then they should 
also apply to the QSO spectra. However, the results in Fig. 2 



II Each set of ThAr lines never included any of the ThAr lines 
from any other set. 
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Figure 2. The top panel shows the results of MOla and W01. The 
middle panel compares the ThAr calibration error on the same 
scale. The lower panel gives an expanded view of the ThAr results. 
The weighted mean of the ThAr sample is (0.0±1.4) X 10~ 7 . ThAr 
data were not available for five of the absorption systems. Clearly, 
this docs not effect our conclusions. 



show that such variations do not produce a systematic effect 
in Aa/a at any detectable level. 

We therefore conclude this section by stating that the 
results in Fig. 2 show unambiguously that wavelength cali- 
bration errors and variations in IP asymmetries are not re- 
sponsible for the observed shifts in a in the QSO absorption 
line results. 



4 DETAILED ANALYSIS: SYSTEMATIC LINE 
BLENDING WITH UNKNOWN SPECIES 

In this section we suppose there is an unidentified transition, 
arising from some species in the same gas cloud as is being 
studied. We further assume that this interloping transition is 
not strong enough and/or sufficiently displaced from one of 
the lines in the analysis, to be detected directly. We describe 
two approaches to explore this: in Section 4.1 we attempt 
to identify candidate blending species and in Section 4.2 we 
describe the results of a test where we remove one transition 
or one species at a time to investigate the impact on Aa/a. 



4.1 Search for potential interlopers 

There are two main parameters which will determine 
whether or not any particular interloper can give rise to a 
significant effect in terms of Aa/a: column density and the 
position in wavelength with respect to the 'host' line of inter- 
est. In order to search for a potential interloper, we explore 
the range in parameter space through numerical simulation, 
generating trial interlopers, blended with a single stronger 



'host' line, and examining the resulting shift in the combined 
centroid in terms of Aa/a. Having established the allowed 
characteristics of a potential interloper, we use photoion- 
ization equilibrium models to search for candidate atomic 
species. We ignore molecular species. 



4-1.1 Characteristics of a candidate interloper 

We restrict the following discussion to blending with Mgn 
A2796, since Aa/a is particularly sensitive to blending with 
this anchor line in the low-z sample. We later generalize 
our discussion to include blending with all other relevant 
transitions. 

Consider an interloper - a spectral line of a blending 
species X - which is separated from the fitted line centroid 
wavelength, Aq, by — AA (AA > 0) in the frame of the cloud. 
Let the separation between the actual Mgn A2796 wave- 
length and the fitted line position be dA (i.e. dA > 0) in the 
frame of the cloud. That is, the line is fit as a single compo- 
nent, despite the presence of the interloper. The interloper 
has a column density N(X). However, as the species X is 
not known a priori, we may only determine the quantity 
N(X)fx where fx is the oscillator strength of the interlop- 
ing transition. 

We generated a synthetic optically thin Mg n A2796 
line and matched this against a composite Mg II A2796 pro- 
file which included an additional weak blended component, 
representing our candidate interloper. The model was fit- 
ted to the 'data' using vpfit. The synthetic spectrum was 
generated using a column density of 4 x 10 12 cm -2 and a b- 
parameter of 5.0 kms -1 as these were representative values 
for the Mgn lines in MOla and W01 (Churchill 1997). 

We restricted the range of 6-parameters for the in- 
terloper using the following physical considerations. If we 
consider all lines to be thermally broadened then the b- 
parameter varies inversely as the square root of the atomic 
mass. Any species lighter than Mg will have a larger b- 
parameter. Since we are considering X to have transitions 
in the ultraviolet, we may assume that it probably has a 
mass greater than or comparable to that of Mg. We there- 
fore place an upper limit of b(X) < 6.0 kms -1 . Also, we did 
not consider any species with atomic number greater than 
~100 and so we adopted a lower limit of b(X) > 0.5 kms -1 . 

By comparing sets of model and synthetic blended lines 
in this way, we can constrain the range of possible values of 
dA, given some apparent Aa/a. 

We present the results of the simulation in Fig. 3 
where we define Fx, Mgn = log 10 [N (X) / N (Mg n)] and Fx = 
log 10 /x- The ripples on the three surface reflect errors 
(5a 0.1 dex in Tx.Mgn) introduced due to the restrictive up- 
per limit placed on b(X). 

Using Fig. 3, we can now specify a certain value of Aa/a 
(i.e. we specify dA) and an interloping transition of a species 
X (i.e. we specify AA) and find the lower bound on the 
column density which that species must have in order to 
have caused the said shift in Aa/a. For example, if we are 
concerned with Aa/a ~ 10~ 5 then dA ~ 50x 10" 4 A and the 
section of Fig. 3 with dA 50 x 10 _4 A contains information 
which restricts the candidate interlopers. 



© 2000 RAS, MNRAS 000, 0-|l| 



8 AT. T. Murphy et al. 




Figure 3. Three space of results from the blending simulation. dX 
directly corresponds to a value of Aq/q if blending is responsible 
for all observed shifts in figure 3 of MOla. The figure is shaded 
according to values of Tx.Mgii- Ripples on the surface reflect er- 
rors due to the restrictive upper bound on b(X). The 'lighting' of 
the figure has been adjusted to maximize their visibility and the 
main errors have been highlighted in black. 

4-1.2 Photoionization Modeling 

Given a species, X, what is an upper limit on the column 
density in any given cloud? To answer this, we modelled the 
absorption clouds with a grid of CLOUDY models (Ferland 
1993). CLOUDY estimates the photoionization conditions in 
absorption clouds by assuming that the cloud can be ap- 
proximated as a series of zones, in a plane parallel geom- 
etry, illuminated by UV radiation from one direction. The 
fact that we are trying to model a (presumably) spherical 
cloud with UV flux incident from all sides is expected to 
introduce abundance errors of order a factor of 2 (Ferland 
1993). 

Photoionization equilibrium is characterized by (a) the 
shape and intensity of the UV spectrum, (b) the number 
density of hydrogen atoms in the cloud, tih, (c) the size 
of the cloud, which we assume is the same for all species 
(the size can be found by specifying rm and iV(Hi), the 
neutral hydrogen column density) and (d) the metallicity, 
Z, of the cloud. It is also useful to define the ionization 
parameter as U = n^/nu, where n 7 is the number density 
of ionizing photons incident on the cloud. Here we follow the 
work in Bergeron & Stasihska (1986), Bechtold et al. (1987), 
Steidel & Sargent (1989) and Bergeron et al. (1994) and 
choose a power law spectrum, /„ oc ^ _3//2 , for /„ the incident 
flux at frequency v. The low and high energy extremes are 
given by /„ oc v 5 ^ 2 (for A„ > 10 (im) and /„ cx v~ 2 (for 
E u > 50 keV). 

All the above parameters (U, nn, JV(Hi) and Z), vary 
from cloud to cloud and are poorly known quantities. We 
therefore use a grid of cloudy models to find an upper limit 
on the column density of a given ionic species using ranges in 
these parameters as follows. For U we adopt the range —4 < 
log 10 U < —3 in line with the measured range of parameters 
for those clouds studied in MOla and W01 (Churchill 1997). 



Similarly, for JV(Hi) we use 17.0 < log 10 JV(Hi) < 19.0. 
For Z we adopt a conservative range, — 2 < Z < (follow- 
ing Morris et al. 1986; Bergeron & Stasihska 1986; Bechtold 
et al. 1987; Steidel & Sargent 1989; Bergeron et al. 1994; 
Churchill 1997). The range of parameters we explore here 
are characteristic of the low redshift sample. For a related 
discussion of DLAs, see Prochaska (1999) and Prochaska & 
Wolfe (2000). 

nu is the most difficult parameter to constrain. To 
obtain an estimate of the lower bound, one may use the 
method, first suggested by Bahcall (1967), whereby one 
places an upper limit on n p , the number density of protons 
in the cloud, using the equation, 

= 3.9 Xl0-V(l+ 0.22^). (5) 

Here, A r (Cn*) is the column density of C II which has 
been collisionally excited to the state Cn*. n e is the elec- 
tron number density. We assume that n e ~ n p and that 
n p ~ n-H, and so, using iV(Cn) = 4.6 x 10 14 cm -2 and 
JV(ClI*) < 10 12 cm -2 from Songaila et al. (1994), we esti- 
mate that log 10 71h < —1.3. Under these simple assumptions, 
this value is consistent with Morris et al. (1986) estimate, 
log 10 riH < 0.0, and we adopt this as a conservative upper 
limit. A lower limit may be obtained by using upper limits on 
the cloud sizes. We used results from Bergeron & Stasihska 
(1986) to place a lower limit of log 10 nu > —3 where we 
have used CLOUDY to estimate that iV(Hl) ~ iV(H)/1000. 
We compare our parameter range estimates above with those 
assumed in Lopez et al. (1999) and find that they are consis- 
tent after considering the differences between the assumed 
radiation field. 

We present the results of the grid of CLOUDY models in 
Table 1 where we only give the value of r™^ gII , the max- 
imum value of rjf,Mgn that occurred in the grid. We also 
present the model responsible for this value of Fx.Mgn- All 
values of F™^ MgII should be treated as having an error of at 
least 0.3 dex. We find that these results are consistent with 
the predictions of model 1 in Lopez et al. (1999) and are 
also consistent with their measured values. 

4-1-3 A Search for Interlopers and Results 

Having obtained estimates of the column density upper lim- 
its (Table 1), we can use the results of Fig. 3 to restrict the 
possible values for AA for any of the species above. Fig. 3 
shows the value of this test, since, given a value for Fx.Mgn, 
we can place very stringent constraints on the allowed range 
in AA. Then, we search tabulated transitions for each species 
to see if any candidate interlopers exist. 

Table 1 and the results expressed in Fig. 3 show that 
any interloper must have a wavelength differing by no more 
than ~ 0.1 A to that of Mgn A2796. We searched published 
energy level tables (Moore 1971; the Vienna Atomic Line 
Database (VALD) - Piskunov et al. 1995 and Kupka et al. 
1999) to find potential lines of those species shown in Ta- 
ble [l] which satisfy the above criteria. We searched only for 
El transitions as all other transitions are associated with a 
large factor of suppression. For example, E2 and Ml tran- 
sition probabilities are suppressed by factors ~ a~ 2 and 
~ or 2 {Z n oi)~ 4 ' respectively (Z n is the nuclear charge). First- 
forbidden El transitions can also be considered since the 
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Table 1. Results of the CLOUDY photo-ionization equilibrium calculation showing the upper limit on Tx.Mgii = l°gio jv(iv%ii) • m °del 
used is specified by the logarithmic parameters, U, iV(Hl), njj and Z. 



Ion 


pmax 
1 X.Mgll 


U 


N(H i) 


"II 


Z 


Ion 


pmax 
1 X,MgII 


U 


AT(Hl) 


"11 


Z 


Fel 


-0.20 


-4 


19 


-3 


-2 


Nil 


-2.39 


-1 


17 


-3 


-2 


Fell 


0.62 


-3 


19 


-3 





Nill 


-0.79 


-3 


19 


-3 





Fc in 


0.84 


-3 


17 


-3 


-2 


Ni III 


-1.07 


-3 


17 


-3 


-2 


Fe iv 


0.99 


-3 


17 


-3 


-2 


Ni iv 


-0.63 


-3 


17 





-2 


Fev 


0.42 


-3 


17 


-3 


-2 


Niv 


-1.05 


-3 


17 


-3 


-2 


Fc vi 


O.lf 


-3 


17 


-3 


-2 


Ni vi 


-0.62 


-3 


17 


-3 


-2 


Fc vii 


-0.71 


-3 


17 


-3 


-2 


Ni vii 


-0.92 


-3 


17 


-3 


-2 


Fe vm 


-2.07 


-3 


17 


-3 


-2 


Ni vm 


-1.79 


-3 


17 


-3 


-2 


Mgi 


-1.50 


-4 


19 


-3 


-2 


Nai 


-2.18 


-1 


17 


-3 





Mgin 


1.16 


-3 


17 


-3 


-2 


Nan 


-0.80 


-3 


17 





-2 


Mg iv 


0.19 


-3 


17 


-3 


-2 


Nam 


-0.23 


-3 


17 


-3 


-2 


Mg v 


-1.24 


-3 


17 


-3 


-2 


Naiv 


-0.74 


-3 


17 


-3 


-2 


Mg vi 


-2.20 


-3 


17 


-3 


-2 


Na v 


-1.82 


-3 


17 


-3 


-2 


All 


-2.02 


-3 


17 


-3 


-2 


Si I 


-1.53 


-4 


17 


-3 





Alii 


-0.17 


-3 


17 


-3 


-2 


Sin 


0.62 


-3 


17 


-3 





Al III 


-0.64 


-3 


17 








Sim 


1.09 


-3 


17 


-3 


-2 


Al iv 


-0.31 


-3 


17 








Si iv 


0.26 


-3 


17 





-2 


Alv 


-1.37 


-3 


17 








Siv 


-0.62 


-3 


17 





-2 


Alvi 


-2.28 


-3 


17 








Si vi 


-1.46 


-3 


17 





-2 


Cai 


-2.03 


-4 


17 


-3 





Ci 


-0.33 


-1 


17 


-3 


-2 


Can 


-1.32 


-4 


17 








Cn 


1.53 


-3 


17 


-3 


-2 


Cam 


-0.01 


-3 


17 


-3 


-2 


Cm 


2.08 


-3 


17 


-3 


-2 


Ca iv 


-0.84 


-3 


17 


-3 


-2 


Civ 


0.85 


-3 


17 


-3 


-2 


Cav 


-2.18 


-3 


17 


-3 


-2 


C v 


0.00 


-3 


17 


-3 


-2 


N i 


0.70 


-4 


19 








Oi 


1.79 


-3 


19 


-3 





Nil 


1.01 


-3 


17 


-3 


-2 


On 


1.88 


-3 


19 


-3 





Nm 


1.57 


-3 


17 


-3 


-2 


Om 


2.47 


-3 


19 


-3 





Niv 


0.48 


-3 


17 


-3 


-2 


Oiv 


1.13 


-3 


19 


-3 





Nv 


-1.03 


-3 


17 


-3 


-2 


Ov 


-0.21 


-3 


19 


-3 





Nvi 


-1.89 


-3 


17 


-3 


-2 


O vi 


-1.70 


-3 


19 


-3 





Nel 


-0.22 


-4 


19 


-3 





An 


-1.24 


-4 


19 


-3 


-2 


Ncn 


0.57 


-3 


19 


-3 





Ar n 


-0.53 


-3 


17 


-3 


-2 


Ne III 


1.68 


-3 


17 


-3 


-2 


Arm 


0.01 


-3 


17 





-2 


Ne IV 


0.68 


-3 


17 


-3 


-2 


Anv 


-0.89 


-3 


17 





-2 


Ne v 


-0.70 


-3 


17 


-3 


-2 


Ar v 


-2.38 


-3 


17 


-3 


-2 



oscillator strengths are typically ~ WZ 4 a 4 times that of a 
normal El transition in a related or similar ion. 

Oscillator strengths, fx, are not generally available for 
the potentially interloping transitions, so we must adopt an 
upper limit by using other (stronger) transitions for the same 
multiplet (for which oscillator strengths are known). 

Given upper limits on fx and Tj^Mgn and a value for 
A A, we can estimate d\ and hence derive a conservative 
upper limit on the interloper's effect on Aa/a. 

Although Fig. 3 relates specifically to Mg II A2796, sim- 
ilar simulations were carried out to search for interlopers 
near all other 'host' absorption lines (i.e. those in table 1 in 
MOla). The model parameters were also extended to DLAs 
(e.g. 19<log 10 iV(Hl)<22). 

We have not been able to identify any candidate inter- 
lopers satisfying all of the criteria above. 

We note that CLOUDY only computes column densities 
for the most abundant species. However, it is possible that 



other low abundance species could mimic a shift in a. There- 
fore, we searched for interloping lines (in Moore 1971 and 
VALD) of ionic species of elements such as Sc, Ti, V, Cr, Mn, 
Co and Zn and found only one possibility: a Cr II transition 
with A = 2026.2686 A which lies 0.13 A to the red of the 
Znn A2026 'host' line. This line has an oscillator strength of 
0.047 according to VALD. We simulated its effect on Aa/a 
for the high- z sample. 

We obtained the largest value of iV(Cr n)/JV(Zn n) di- 
rectly from the observed QSO absorption line data. We then 
generated a synthetic, single velocity component spectrum 
containing Zn II A2026 and the Cr n interloper. We also gen- 
erated several other lines seen in the high-z sample. We pro- 
duce synthetic spectra with several different values of Aa/a 
and then used vpfit to measure Aa/a. The result was that 
the measured value of Aa/a was always consistent with the 
nominal value, that is, we could detect no effect due to the 
interloping Cm transition. 
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In summary, we have found no interloping species which 
could have caused a significant shift in Aa/a. We do note, 
however, that the ionic energy level data in Moore (1971) 
and VALD are incomplete for many highly ionized species 
and so we cannot rule out the possibility of blending with 
undiscovered transitions. Also, as noted previously, this 
analysis only takes into account blends with atomic, not 
molecular species. Firstly, typical oscillator strengths for 
molecular transitions are much smaller than for resonance 
atomic transitions. Also, the typical column densities of the 
low abundance elements in Table |l| suggest that very few 
molecules will exist in the QSO absorption clouds. There- 
fore, we consider the possible effect of any molecular blend- 
ing to be negligible. 

4.2 Line removal as a test for systematic blends 

If a particular host transition (e.g. Mgn A2796) were sys- 
tematically blended then removing that transition from the 
analysis in MOla and W01 should result in a different value 
of Aa/a. Line removal is also a rough test for isotopic ra- 
tio or h yperfine structure effects as introduced in Sections 
2.3-2.5. We therefore remove a particular transition in an 



absorption complex only if doing so still allows us to obtain 
meaningful constraints on Aa/a. For example, we do not 
remove the Mgn A2796 line from a fit to a Mgll/Fell sys- 
tem if the Mg II A2803 line is not present. We also removed 
entire species where possible. A small difficulty occurs with 
Cm and Znll A2062, whose profiles are sometimes blended 
together (i.e. when the velocity structure is broad). If the 
profiles were blended then we simply removed both transi- 
tions simultaneously, only removing them separately when 
the lines were completely separated. We obtain Aa/a with 
the same procedure described in MOla. 

Removing a transition from the fit may result in a 
slightly modified estimate for the velocity structure, revis- 
ing the estimate for Aa/a for that system. The question 
is, how robust is the estimate for Aa/a, averaged over the 
entire sample, to this line removal process? 

We present our results in Fig. 4. For each transition re- 
moved, we compare the weighted mean of Aa/a (and the 
associated la error) before and after the line removal. It 
must be stressed that both the value and the errors are not 
independent and so interpreting Fig. 4 is difficult. However, 
it seems that there are no extreme deviations after line re- 
moval. Thus, we might conclude from such a diagram that 
systematic blending has not occurred and that we have no 
reason to suspect that one or more of the laboratory wave- 
lengths is incorrect. 

Concentrating on the heavier ions, we might also con- 
clude that there is no effect on Aa/a due to possible iso- 
topic ratio or hyperfine structure effects. We cannot draw 
similar conclusions for Mg since we cannot remove all Mg 
lines simultaneously (this would prevent us from obtaining 
meaningful constraints on Aa/a for the low-z sample). We 
also see that any lack of thermal equilibrium in the absorp- 
tion clouds, leading to un equa l population of the hyperfine 
levels of Aim (see Section 2.5), has had no significant effect 



on Aa/a. The large error bars on the Aim points reflect 
the fact that those lines are only present in 6 of our high- 
z systems. 

Qualitatively, this line removal test suggests that there 
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Figure 4. Comparison of the weighted mean of Aa/a before 
(dotted error bar) and after (solid error bar) line removal. The 
transitions removed are listed on the left together with the num- 
ber of systems, n, for which removal of that transition was possi- 
ble. It is stressed that comparing the values and error bars before 
and after line removal is difficult since these quantities are not 
independent. Note that there may be some confusion due to the 
occasional blending of the Cm and Znll A2062 lines: '(Cm + 
Znll) A2062' refers to cases where both transitions had to be re- 
moved simultaneously; 'Cm + Znll A2062' refers to cases when 
all Cr II transitions were removed along with the blended Zn II 
line; a similar definition applies to 'Znll + Cm A2062'; 'Cm' 
then refers only to the removal of all Cm transitions in cases 
where there was no blending with the Znll line. Only one simi- 
lar case occurred for removal of all Znll lines and so we do not 
present this result. 



is no one single 'host' transition or species which is system- 
atically blended and which has therefore affected our results 
in any significant way. The results also suggest that isotopic 
ratio or hyperfine structure effects have not significantly ef- 
fected the heavier ions in our analysis. This point is more 
fully explored in Section | and §. 



5 DETAILED ANALYSIS: ATMOSPHERIC 
DISPERSION EFFECTS 

The problems that arise due to atmospheric dispersion ef- 
fects can be understood with reference to Fig. 5. If we ob- 
serve a point source with the spectrograph slit at some angle 

6 to the vertical (as projected on the sky) then the object is 
dispersed with some component along the spectral direction 
of the slit. The optical design of the Keck/HIRES combina- 
tion at the time at which a large fraction of our data was 
taken was such that 6 — £ for £ the zenith angle of the ob- 
ject (Tom Bida, Ste ve Vo gt, personal communication). As 
we noted in Section 2.12, this leads to two effects on the 
wavelength scale: 
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To zenith 



Seeing discs 
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Green 
Red 




- Slit jaws 
Horizon line 



Figure 5. Schematic diagram showing the effects of atmospheric 
dispersion. The HIRES slit is projected onto the sky so that we see 
the seeing discs of blue, red and green light from a point source 
dispersed with some component along the spectral direction of 
the slit. The spectral direction is indicated by A and the spatial 
direction by X. Note that the refraction is worse in the blue. 



(i) The first effect is a stretching of the spectrum (relative 
to the ThAr calibration frames) due to the angular separa- 
tion of different wavelengths entering the slit. Consider two 
wavelengths, Ai and A2 (A2 > Ai), falling across the slit as 
in Fig. 5. If we were to measure the spectral separation be- 
tween these two wavelengths on the CCD, AA', we would 
find that 

aAtp sin 6 



AA' » A 2 - Ai + 



(6) 



where a is the CCD pixel size in An gstroms, S is the 
projected slit width in arc seconds per pixel (for HIRES, 
5 — 0'.'287 per pixel) and Alp is the angular separation (in 
arc seconds) of Ai and A2 at the slit. Atp is a function of the 
atmospheric conditions along the line of sight to the object 
and can be approximated using the refractive index of air 
at the observer and the zenith distance of the object (e.g. 
Filippenko 1982). Note that equation ^ assumes that the 
seeing profile is not truncated at the slit edges. If tracking 
errors or seeing effects cause truncation, equation [] provides 
an upper limit to the measured AA'. 

(ii) If truncation of the seeing profile does occur, it may 
be asymmetric (as the diagram illustrates). In the case of 
HIRES, the optical design is such that a blue spectral feature 
will have its red wing truncated and a red spectral feature 
will be truncated towards the blue. This is the case shown 
in Fig. 5. Thus, an asymmetry is introduced into the IP and 
this will be wavelength dependent. The extent of the asym- 
metry will depend on Aip (which depends most strongly on 

0- 

Note that when we centroid spectral features to find 
their separation, both effect (i) (shifting) and effect (ii) (in- 
strumental profile (IP) distortion) effectively 'stretch' the 
spectrum. For the low-z sample, this should lead to sys- 
tematically positive values of Aa/a. That is, atmospheric 
dispersion effects cannot explain the average negative value 




30 40 50 

Zenith angle (degrees) 

Figure 6. The lower panel shows the distribution of absorption 
systems over average zenith distance. The upper panel shows the 
QSO data as a function of zenith distance. We do not see any ob- 
vious correlation. When we remove points for which there should 
be no atmospheric dispersion effect, the Spearman rank corre- 
lation coefficient is 0.18 with an associated probability of 0.41. 
The coefficients for the low (filled circles) and high (crosses) red- 
shift samples are 0.16 and -0.09 with probabilities 0.43 and 0.71 
respectively. 



of Aa/a seen in the QSO data. Removing these effects from 
our data would make the results more significant. Although 
there is no obvious correlation between Aa/a and £ as would 
be expected in this scenario (see Fig. 6), the potential effect 
is quantified in our analysis below. 

For all objects which were observed without the use of 
an image rotator, we have used the observational parameters 
reported in the image header files to find average values of £ 
(see Fig. 6) , temperature, pressure and relative humidity for 
each object. This allows us to correct for the stretching of 
the spectrum implied by equation ^| We have re-calculated 
Aa/a in the same manner as described in MOla and we plot 
the results in Fig. 7. For several QSOs, only a single image 
header (i.e. for one of the QSO exposures out of a series) was 
available. In these cases we have added (in quadrature) an 
error of 0.9 x 10~ 5 . This is an average value determined by 
calculating Aa/a for several values of the zenith distance for 
those objects concerned and corresponds to an uncertainty 
in the average £ of about ±20° . 

We have also tested the effect of the IP distortions (ef- 
fect (ii) above) on Aa/a. We generated high SNR, single 
component absorption spectra (with Aa/a — 0) for sev- 
eral redshifts. We included different combinations of tran- 
sitions in the different spectra so as to make a representa- 
tive estimate for both samples of QSO data (particularly 
the high-z sample which contains many combinations of dif- 
ferent transitions). A wavelength dependent PSF was con- 
structed using the following procedure. We assume that the 
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Table 2. Comparison of the weighted means, (Aa/a) w , with 
the results of MOla and W01, after correction for atmospheric 
dispersion effects. 



Redshift 



1 



10 



Sample 


<Aa/a> w /10- 5 


(Aa/a)^ 01a /10- 5 


Low- 2 


-1.52 ± 0.23 


-0.70 ±0.23 


High-z 


-0.79 ± 0.26 


-0.76 ±0.28 


Total 


-1.19 ± 0.17 


-0.72 ±0.18 



seeing (taken as 0'.'8) and tracking error (taken as 0V25) gen- 
erate Gaussian profiles across the slit and that a wavelength 
Aq falls at its centre. Other wavelengths are refracted to dif- 
ferent parts of the slit depending on the zenith distance. To 
compute this, typical Keck atmospheric conditions (temper- 
ature, pressure, relative humidity) were used. 

We truncate the profile by multiplying by a top-hat 
function (whose width is equal to that of the slit) and 
then convolve the result with a Gaussian IP of width a — 
2.2 kms -1 (Valenti et al. 1995). This final PSF is then con- 
volved with the synthetic absorption spectra. The value of 
Ao is ~ 5500 A (i.e. the TV acquisition system is centred at 
about this wavelength). Using this value for Ao and £ = 45° 
we find Aa/a to be quite insensitive to the IP distortions, 
particularly in the high-z regime. The effect on Aa/a may 
be as large as 0.05 x 10 for the low-z regime but round-off 
errors in vpfit become important at this level of precision. 
To exaggerate the effect, we also chose extreme values of Ao 
to be 4000 and 7000A. The effect on Aa/a in these cases 
was still negligible in the high-z regime but increased to the 
level of Aa/a ~ 0.12 x 10 -5 for the low-z spectra. This lat- 
ter value corresponds to ~l/6 of the correction due to the 
shifting effect. The sign of the above corrections to Aa/a 
are such that applying them to the QSO data make the ob- 
served average Aa/a more negative. 

In summary, we find Aa/a to be most sensitive to the 
shifting effect described by equation ^ and reasonably in- 
sensitive to the IP distortion effect. Thus, we do not apply 
a correction for the latter in the results shown in Fig. 7. 
The results show the values of Aa/a once the shifting ef- 
fect has been removed from the data. As expected, we see a 
significant decrease in Aa/a in the low-z sample. The high- 
z sample is insensitive to this systematic effect, as can be 
expected due to the many different transitions used with qi 
coefficients of different signs. We summarise the results in 
Table ^| and compare them with those found in MOla and 
W01. We stress that the corrections above (equation ^) are 
upper limits. Seeing and tracking errors will diminish the 
effects substantially. 



6 DETAILED ANALYSIS: DIFFERENTIAL 
ISOTOPIC SATURATION 

In this section we quantify the effect of differential isotopic 
saturation on our measured values of Aa/a. As discussed 
in Section 2.3, this comes about when the strongest isotopic 



component of an absorption line saturates. The weaker iso- 
topic components may still be on the linear part of their 
curves of growth. Using a single weighted wavelength to 
represent the composite profile position (i.e. the profile cen- 



c 
< 
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-1 
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-3 
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Figure 7. Results after removal of atmospheric dispersion effects. 
The top panel shows our raw results for Aa/a as a function of 
look-back time in a flat, A cosmology. The redshift scale is also 
provided for comparison. The lower panel shows an arbitrary bin- 
ning of the results which emphasizes the susceptibility of the low- 
z sample to this systematic error (the dotted lines are the results 
of MOla and W01, slightly shifted for clarity). Note, however, that 
the correction made here is an extreme and will be diminished by 
seeing and tracking effects. 



troid) can therefore, in principle, result in systematic errors 
in Aa/a. 

The various isotopic abundances of concern are shown 
in Table ^. If only composite wavelengths were used then, if 
there exists a saturated dominant isotope in the QSO data, 
this may result in systematic errors in Aa/a. However, com- 
posite wavelengths were only used for the 4 heaviest atoms 
(i.e. Cr, Fe, Ni, Zn) and, of these, only some Fe transitions 
appear to be saturated in some of our QSO spectra. The Cr, 
Ni and Zn lines were always very weak (typically absorbing 
only ~20% of the continuum) and so using composite tran- 
sition wavelengths is justified. 

It is therefore left to quantify the effect of differential 
isotopic saturation for the Fe lines in our spectra. Ideally, if 
we knew the isotopic wavelength separations for the various 
Fe II transitions then we could do this by simply fitting our 
spectra with this structure to find Aa/a. The difference 
between the weighted mean Aa/a and that found in MOla 
and W01 would indicate the size of the effect. In reality, 
we do not know the isotopic separations. We can, however, 
place an upper limit on this effect since we know (or can 
estimate) the isotopic separations of the Mg and Si anchor 
lines. If we fit our spectra with only the composite Mg and Si 
wavelengths then we will obtain maximal changes to Aa/a 
since (i) these atoms have strong lines which saturate more 
often than the Fell lines and (ii) they will have the largest 
isotopic separations (i.e. largest mass shift). 

We have carried out the above re-calculation of Aa/a 
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Table 3 . The percentage terrestrial abundances of the atoms used 
in our analysis (Rosman & Taylor 1998). The neutron number is 
defined relative to that of the isotope with the largest abundance 
and is denoted as AA (negative values representing lighter iso- 
topes). 



Atom 


-2 





+1 


AA 

+2 


+3 


+4 


+6 


Mg 




78.99 


10.00 


11.01 








Al 




100.0 












Si 




92.23 


4.68 


3.09 








Cr 


4.345 


83.789 


9.581 


2.365 








Fe 


5.845 


91.754 


2.119 


0.282 








Ni 




68.08 




26.22 


1.14 


3.64 


0.93 


Zn 




48.6 




27.9 


4.1 


18.8 


0.6 



using the composite wavelengths for the Mgl, Mgn and 
Sill lines. In the low-z sample, the fact that the Mg lines 
have their weaker isotopes at shorter wavelengths means 
that Aa/a should be systematically more negative if the 
composite wavelengths are used. The situation is again 
more complicated in the high-z sample. We find only very 
small differences between the results and those of MOla 
and W01, the weighted means for the low and high-z sam- 
ple being Aa/a = (-0.75 ± 0.23) x 10~ 5 and Aa/a = 
(-0.74 ± 0.28) x 10 -5 respectively (compared to Aa/a = 
(-0.70 ± 0.23) x 10~ 5 and Aa/a = (-0.76 ± 0.28) x 10~ 5 
from MOla and W01). 

Can we understand the magnitude of these corrections? 
Imagine a single velocity component and consider finding 
Aa/a from, say, the Mgn A2796 and Fell A2344 lines only. 
One might estimate the maximum effect of differential iso- 
topic saturation by assuming that when all isotopes of Mg II 
are saturated, the line centroid would lie at the unweighted 
mean isotopic wavelength. This would result in a correction 
to Aa/a of ~ — 1 x 10 -5 (i.e. the observed value would 
become more negative when the correction is applied). 

However, several factors will reduce this correction in re- 
ality. Firstly, once a line begins to saturate, then because the 
intrinsic line width (i.e. the ^-parameter) is far greater than 
the isotopic separations, the weaker isotopes are swamped 
by the dominant one. This reduces the effect on Aa/a. Sec- 
ondly, if the velocity structure is more complex (i.e. if it 
contains more than one velocity component) the unsatu- 
rated components, which do not suffer from the differen- 
tial isotopic saturation problem, will provide much stronger 
constraints on Aa/a. Finally, as more unsaturated transi- 
tions (e.g. Fe n) are incorporated into the fit of a particular 
QSO absorption system, the effect is further reduced since 
some constraints on Aa/a come from shifts between the 
Fe II lines. 

It is therefore clear that differential saturation of the 
isotopic components of the various Fe II lines cannot have 
lead to any significant effect on Aa/a. Our results also in- 
dicate that the weighted mean Aa/a for the high-z sample 
will not be very sensitive to errors in our estimate of the Si II 
isotopic separations. 



Table 4. Wavelengths of the dominant isotopic components of the 
Mg and Si transitions. Those of Mg are derived from experimental 
data in Pickering et al. (1998) while those of Si are calculated in 
MOla on the basis of a scaling by the mass shift from the basic 
Mg isotopic structure. 



Transition 


rn 




Ao 


Mgl A2853 


24 


35051.271 


2852.9636 


Mgn A2803 


24 


35669.286 


2803.5324 


Mgn A2796 


24 


35760.835 


2796.3553 


Sin A1808 


28 


55309.330 


1808.0132 


Sin A1526 


28 


65500.442 


1526.7073 



7 DETAILED ANALYSIS: ISOTOPIC 
ABUNDANCE VARIATION 

The analysis we have carried out so far assumes terres- 
trial isotopic abundances. However, it is quite possible, even 
likely, that these high redshift gas clouds have rather differ- 
ent abundances. As we noted in Section 2.4, a substantial 



systematic shift in a might be mimicked and will be partic- 
ularly pronounced in the low-z sample since the Mg isotopic 
separations are the largest. No observations of high redshift 
isotopic abundances have been made so we have no a pri- 
ori information on the QSO absorption systems themselves. 
However, observations of Mg (Gay & Lambert 2000) and 
theoretical estimates for Si (Timmes & Clayton 1996) in 
stars clearly show a decrease in isotopic abundances with 
decreasing metallicity. For example, at relative metal abun- 
dances [Fe/H] 1, 25 Mg/ 24 Mg ~ 26 Mg/ 24 Mg w 0.1, i.e. 

about 20% below terrestrial values (see Table B). Theoreti- 
cal calculations (Timmes, Woosley & Weaver 1995) suggest 
even larger decreases. 

Estimates have also been made of the metal abundances 
in a subset of the QSO absorption systems used in our 
analysis. The DLAs have a mean metallicity of Z ~ —1.5 
(Prochaska & Wolfe 2000) and the Mg/Fe systems span the 
range Z = -2.5 to 0.0 (Churchill et al. 1999). 

Therefore, we explored the effect of isotopic abundance 
evolution on Aa/a for the entire range in isotopic abun- 
dance from zero to terrestrial values (i.e. in the case of Mg, 
the limiting case is 25 Mg/ 24 Mg = 26 Mg/ 24 Mg = 0). This 
allows us to derive a secure upper limit on the magnitude of 
the effect on Aa/a. 

We have re-computed Aa/a for the entire sample af- 
ter removing the following isotopes: 25 Mg, 26 Mg, 29 Si and 
30 Si. The wavelengths for the dominant isotopic compo- 
nents, 24 Mg and 28 Si, for the various transitions of interest, 
are given in table 3 of MOla and are summarized in Table 

The results are presented in Table H and Fig. 8. These 
show, as expected, that the low-z sample is more sensitive 
to this systematic error than the high-z sample. We stress 
that the signs of the above corrections are such that ap- 
plying them to the QSO data make the observed average 
Aa/a more negative. Since we have no information on the 
isotopic separations for the other atoms of interest, we can- 
not perform a similar analysis for all species. However, the 
corrections due to variation in the isotopic abundances of Mg 
and Si should be dominant (as we discuss is Section [2,4|) so 
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Table 5. Comparison of the weighted means, (Aa/a) w , with the 
results of MOla and W01, after maximal correction for possible 
variation in isotopic ratios 



Sample 


<Aa/a> w /10- 5 


(Aa/a)^ 01a /10- 5 


Low z 


-1.01 ± 0.23 


-0.70 ±0.23 


High z 


-0.87 ±0.26 


-0.76 ±0.28 


Total 


-0.96 ± 0.17 


-0.72 ±0.18 



Redshift 
1 2 3 




Fractional look— back time 



Figure 8. Results after removing the weaker isotopic compo- 
nents of Mg and Si. The upper panel shows our raw results and 
the lower panel shows an arbitrary binning which emphasizes the 
susceptibility of the low-z sample to this systematic error (the 
dotted lines are the results of MOla and W01, slightly shifted for 
clarity). Note, however, that the correction made here is extreme. 

the corrections shown in Table ^ and Fig. 8 are conservative 
upper limits. 



8 A GENERAL TEST FOR SYSTEMATIC 
ERRORS 

Although we have already convincingly ruled out any signif- 
icant wavelength-scale distortion using ThAr spectra (Sec- 
tion ^), the QSO and ThAr light does not follow precisely 
the same light path. Indeed, for these reasons we had to 
consider some of the previous effects such as atmospheric 
dispersion. Thus, it is conceivable that some residual long- 
scale wavelength distortion or other simple systematic error 
remains in the QSO spectra. 

The dependence of Aa/a on the gi coefficients is simple 
for the low-z sample but considerably more complicated for 
the high z sample. This can be seen by inspecting table 1 in 
MOla. If systematic errors are responsible for the non-zero 
Aa/a we measure, it is therefore surprising that Aa/a is 
so consistent between the 2 samples. A priori, it would be 



Table 6. Results of a general test for systematic errors. We re- 
moved the three different types of transition from a sub-sample of 
the 21 high-z systems. The new Aa/a is given for each case. These 
should be compared to the value obtained when only transitions 
with mediocre q\ coefficients were removed. Note that neither the 
values or the la errors are independent of each other. 



Sub-sample removed 


(Aa/ a ) w /10- 5 


None 


-0.90 ±0.34 


Mediocre-shifters 


-1.29 ±0.39 


Anchors 


-1.46 ±0.44 


Negative-shifters 


-1.34 ±0.65 


Positive-shifters 


-1.39 ± 1.01 



even more surprising if we were to take subsets of the QSO 
data, grouped according to the sign and magnitude of gi, 
and found a consistent value for Aa/a. 

We apply this idea to the high-z sample only, since it 
contains transitions with a wide range of gi values. Of the 21 
absorbers in this sample, we can find a subset of 13 systems 
that include at least one anchor line (transitions with small 
gi), at least one positive-shifter (gi > 700 cm - ) and at 
least one negative-shifter (gi < —700 cm -1 ). If some kind 
of general long-scale distortion of the wavelength scale was 
present, we might expect to find a change of sign in the 
average Aa/a when comparing, for example, the negative- 
shifters or positive-shifters against the anchors (see equation 

0). 

From these systems we remove all transitions with gi co- 
efficients of 'mediocre' magnitude (300 < |gi| < 700 cm -1 ) 
so as to clearly delineate the three types of transition (an- 
chor, positive and negative-shifter). We then calculate Aa/a 
for this sample. We then remove each type of transition sep- 
arately and calculate Aa/a again. There is a small compli- 
cation due to the presence of the occasional blend between 
Cr II and Zn II A2062 since these lines shift in opposite direc- 
tions. In cases where these two lines were blended together, 
we removed both of them from all fits since they are of dif- 
ferent 'type' and cannot be removed individually. 

We present our results in Table |(| Again, precise, quan- 
titative comparison of the values of Aa/a obtained with the 
various types of transitions removed is difficult since neither 
the values or the la error are independent. Also, when we re- 
move the positive shifters, the uncertainly in Aa/a is quite 
large. However, we do not find any evidence to suggest that 
Aa/a is somehow dependent on the sign of gi from this 
test. That is, the high-z results seem to be robust against 
the presence of simple systematic errors. Indeed, this is cor- 
roborated by the results of removing atmospheric dispersion 
effects (Table | and Fig. 7). 



9 CONCLUSIONS 

The most important conclusion of the work described here 
is that after a detailed search for instrumental and astro- 
physical effects, we have found no systematic effect which 
can easily mimic a negative Aa/a. 

Of those effects we have identified and explored, two sys- 
tematic effects have been found to be potentially significant: 
atmospheric dispersion and isotopic evolution. If we were 
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to remove these effects from the data, the resulting Aa/a 
would actually be more negative than the uncorrected effect 
(see Figs. 7 and 8). These corrections are not applied to the 
data in MOla and WOf because (i) our estimates provide 
upper limits on the magnitude of the effects, and (ii) we pre- 
fer to be as conservative as possible as to the significance of 
the deviation of Aa/a from zero. 

Tighter constraints on a are possible for realistic SNR 
from current telescopes/spectrographs. The analyses above 
indicate that another order of magnitude precision can be 
gained using our methods before wavelength calibration and 
laboratory wavelength errors become significant. We hope 
that MOfa, W01 and the present work encourage further 
well calibrated, higher SNR observations of QSOs and that 
an independent check on our results can be made. 
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